Application of the scattering rate sum-rule to the interplane optical conductivity of 
high temperature superconductors: pseudogap and bi-layer effects 
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We use a recently proposed model of the interplane conductivity of high temperature super- 
conductors to investigate the "scattering rate sum-rule" introduced by Basov and co-workers. We 
present a new derivation of the sum-rule. The quantal and thermal fluctuations of the order param- 
eter which have been argued to produce the observed pseudogap behavior are shown to increase the 
total integrated "scattering rate" but may either increase or decrease the "quasiparticle" contribu- 
tion from frequencies greater than twice the superconducting gap. 
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One of the crucial current questions in the field of high 
temperature superconductivity is the physics of the pseu- 
dogap regime observed in underdoped systems 0]. This 
regime is characterized by the opening of a gap in the 
electronic excitation spectrum. The origin of the gap 
is the subject of debate, with some authors attributing 
it to a density wave of one sort or another Q and oth- 
ers attributing it to the superconducting gap, with long 
ranged superconducting order prevented by loss of phase 
stiffness A second, possibly related, issue is the de- 

gree to which the superconducting ground state differs 
from the familiar BCS mean field state. In particular, 
the possibility of strong superconducting phase fluctua- 
tions signaling proximity to a nontrivial insulating phase 
has attracted considerable interest 

Measurements of the optical conductivity are an im- 
portant source of information about the pseudogap 
regime. For example, the absorptive part of the inter- 
plane [a- axis) conductivity shows a dramatic suppres- 
sion at low frequencies as the pseudogap is formed. The 
detailed frequency and temperature dependence was re- 
cently argued to support a superconducting origin of the 
pseudogap Very recently, Basov and co-workers 

have introduced a new way of interpreting optical con- 
ductivities 0], based on a "scattering rate sum rule" in- 
volving the real part of the inverse of the complex con- 
ductivity which is defined by 



Re a 



Re <t2 + Im 0-2 



(1) 



In this paper we provide a derivation of the sum rule 
which is convenient for our analysis and analyze the the- 
ory of refs [^y^ for single-layer high-Tc superconductors 
such as La2-xSrxCuOi and it's recent generalization | [T^ ] 
to hi-layer materials such as Y Ba2Cu30e+x in terms of 
the new scattering rate sum rule. We show that the so- 
defined "scattering rate" behaves in a nonintuitive man- 
ner, increasing in the pseudogap regime. This behavior 
may be a valuable probe of the quantal and thermal phase 
fluctuations of the superconducting order parameter. 



We first present an alternative derivation of the scat- 
tering rate sum-rule. We suppose that at large frequency 
and for any temperature, Recr^^(aj oo, T) — > a^^ with 
CTo a constant independent of uj and T. In the famil- 
iar Drude model, <Tq^ is proportional to the scattering 
rate, while in realistic models electrons decouple from 
the underlying crystal at sufficiently high frequencies, so 
ctq"^ — 0. This assumption implies that the inverse con- 
ductivity a^^{uj,T) which is a causal response function, 
obeys a "subtracted" Kramers-Kronig relation. 



Ima [lu) — — 
J-oo 7^ 



dx Re (T ^ (x) — (Tq 



UJ ~ X 



(2) 



We use the w ^ oo limit of Eq g to obtain an equa- 
tion for the changes in the inverse conductivity as the 
temperature is changed from Ti to T2, 

lim a; (lmcr^"'^(ct;, Ti) — Inio'^"'^(cij, r2)) 

f°° 2dx 

= / Re{a~\x,Ti)-a-\x,T2) . (3) 

Jo ^ 

In many circumstances, including the c-axis conductivity 
of high temperature superconductors, the integral on the 
right hand side converges rapidly with frequency. 

We now apply this sum rule to a theoretical model re- 
cently proposed |5|j^ for the pseudogap behavior of the c- 
axis conductivity of underdoped high temperature super- 
conductors. In these papers it is argued that the weak- 
ness of the interplane coupling allows one to treat the 
problem by perturbation theory in the interplane cou- 
pling, leading to an expression for the c-axis conduc- 
tivity involving a convolution of normal and anomalous 
in-plane Green functions. The model has three param- 
eters: an in-plane scattering rate F, a superconducting 
gap A, and a Debye Waller factor a characterizing the 
strength of in-plane thermal and quantal fiuctuations of 
the phase of the superconducting order parameter. The 
observed rf-wave structure of the superconducting gap 
Q and the variation around the Fermi surface of F |^ 
may be neglected because the interplane coupling is only 
appreciable in momentum regions in which both F and 
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A are large ||T^,|Tl[. Photoemission ||] and c-axis opti- 
cal data show that in the relevant region of the Fermi 
surface, F is very large and has negligible temperature 
dependence. We take the F ^ oo limit in what follows 
so that in the model of refs temperature enters via 
the values of A and a. At high temperatures, above the 
pseudogap regime, A = a = 0. The pseudogap regime is 
characterized by a A 7^ but a = (superconductivity 
destroyed by phase fluctuations). In the superconduct- 
ing state Of 7^ 0, and an important question is whether as 
T ^ 0, a — > 1 (as in BCS mean field theory) or a < 1, 
signaling quantum fluctuations. The results of refs [^] 
lead to 

T) ^ao[l + iAiT {a - 1) /2uj + 0{l/uj^)] . 

(4) 

Combining Eqs ^ and ^ yields our final expression 
f°° 2dx 

/(a) = / Re(cr-i(a;,A,a)-cr-^(x,A = 0,a)) 

Jo 

_ A7r(l - a) 



Thus the opening of the pseudogap (A > 0, a = 0) cor- 
responds to an increase in the area under the Rea^^ 
curve, while the subsequent entry into the superconduct- 
ing state (a > 0) leads to a decrease. 

The sum rule on Rea"^ is sometimes interpreted 
as a "scattering rate sum rule", because it is conven- 
tional to write a in the 'modified Drude form' a{Lu) = 
LUp/ {—iujX{uj) + in which the real part of a^^ is 

interpreted as a scattering rate. Because the opening of 
the pseudogap, which corresponds to a decrease in the 
density of states, leads to an increase in the integrated 
area of the Recr~^ curve, care is required in using this 
interpretation of scattering rate sum rule data. 

We now apply the results of ||J^ to study in more de- 
tail the behavior of Recr^^ in the pseudogap regime of 
a single-layer superconductor. Fig 1 plots the contribu- 
tions to = ReCT^^(w, A, a) for the normal state and 
for superconducting state at u > 2A for several values of 
a. For a < 1, A 7^ 0, there is an additional delta function 
contribution to arising from a zero-crossing of Imcr 
at < 2 A (not shown in Fig 1). The dashed flat curve 
(labeled by n) corresponds to the normal (A = 0) state, 
while the other curves correspond to the superconduct- 
ing state for different values of a. One sees that the form 
and magnitude of t"-'^ depend sensitively on the strength 
of the fluctuation parameter. For some a-values (0.7 for 
example) the area under for to > 2A is manifestly 
greater than for the normal state, whereas for small a 
(a = 0.1 for example) it is manifestly less. 
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Fig. 1 Optically defined scattering rate (Eq |i|) for single- 
layer superconductor normalized to a; — > 00 limit for normal state 
(dashed line) and for superconducting state with different values 
of fluctuation parameters a = 0.1, 0.7 and 1. For superconducting 
state only the contribution for u> > 2A is shown. 

To further characterize the structure in r^^ it is useful 
to define 

r—Rea-\x,A)- r —Rea-\x,A = 0) (6) 

J2A Jo 

which excludes from Eq ^ the contribution coming from 
any w < 2A delta functions. This quantity is shown 
in Fig 2 along with the normalized total area difference 
f{a). One sees that for a near 1 the increase in 
comes mainly from the lu > 2A region, whereas as a — > 
there is a decrease in the w > 2 A contribution to , but 
the increased value of the 'pole' contribution at w < 2A 
leads to a net increase in the sum. 

1.5 I ' , ' , ' , ' , ' 1 



1 - 




0.2 0.4 0.6 0.8 1 

1-a 

Fig. 2 Difference in the integrated area of (Eq |l|) between 
the gapped (superconducting or pseudogap) state and ungapped 
(A = 0) state as function of fluctuation parameter a. Open circles: 
quasiparticle (i^ > 2 A) contribution only (Eq|6|). Dotted line : to- 
tal contribution. Difference f(a) — /ij>2A(f*) is greater than zero 
for all a < 1, but is too small to be seen for a > 0.7. 

To elucidate this behavior we show in Fig 3, the real 
and imaginary parts of the conductivity for a = 1 (fully 
phase coherent or BCS), a = 0.6 and 0.1 superconductor. 
For a — I (Fig 3a) the imaginary part diverges as 1/lu 
as w — > and with increasing uj monotonically goes to 
zero (as (jJ~^ at high frequencies). The real part increases 
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gradually from zero at the gap edge due to the type-II 
coherence factors. The in this case is zero below the 
gap and is a smooth function above the gap as shown in 
Fig 1. The weak maximum in is close to the position 
of a minimum in Hea^+ Imcr^. 
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Fig. 3 Real and imaginary parts of tlie conductivity for 
a = 1,0.6 and 0.1 superconducting state of single-layer supercon- 
ductor calculated as described in refs ■ 

Fig 3b shows Re and Im a for a typical 'intermedi- 
ate' fluctuation parameter a = 0.6. The crucial effects 
of reducing a from unity are seen to be (1) elimination 
of the type-II coherence factor effect at the gap edge (so 
within this model the Recr acquires a step at the gap 
edge) and (2) a corresponding negative divergence of Imcr 
at a; = 2A, leading to two zero-crossings near 2A. The 

> 2A zero crossing combined with the smallness of Re 
a gives rise to a sharp peak in t^^ (appearing to be at the 
gap edge on the scale of the plot in Fig 1) and the uj < 2 A 
zero crossing leads in this model to a weak delta-function 
contribution to as Retr is zero in that region. At a 
higher frequency the Im part crosses over to the negative 
side but r^^ continues to decrease monotonically. Fig 3c 
shows Re and Im a for strong fluctuations, a = 0.1. The 
increase in fluctuations substantially weakens the l/w di- 



vergence in Imcr and strengthens the negative divergence 
at the gap edge, eliminating the lu > 2A zero crossing 
and moving the to < 2A zero crossing to lower frequen- 
cies, substantially strengthening the delta function at low 
frequencies. In the pseudogap regime (a = 0, not shown) 
the oj < 2A zero crossing (and the associate delta func- 
tion in ) moves to = 0. 

We will now outline the new features that are intro- 
duced in the analysis of Rea"^ when we consider the 
generalization of to the bi-layer case. A detailed 
discussion of the c-axis conductivity for the bi-layer case 
is given elsewhere The essential point is that in a 

bi-layer system a uniform applied field may lead to a non- 
uniform charge distribution within a unit cell which in 
turn produces electrochemical potential differences which 
affect charge flows. These "local-field effects" give rise for 
example to the bilayer plasmon feature observed in many 
bi-layer materials such as YBa2Cu30Q+x fl^ - 

We consider here a bi-layer structure with a unit cell 
consisting of two identical superconducting planes sep- 
arated by a distance di, and with a plane in one cell 
separated from the nearest plane in another unit cell by 
distance c?2- The approach used in ||^,^ yields 



^bilayer 



(Tl(T2 - iuj{cridl -\- (T2d2)/C 
(TlC?2 + 0'2di — iijj/C 



(7) 



with di^2 = di^2/{di + ^2) and cri,2 the conductivities ap- 
propriate to a pair of planes separated by distance di^2 
and C a blockade parameter expressing the electrochem- 
ical potential difference created by a charge imbalance 
within the unit cell. We assume the two conductivities 
to be the same up to a scale factor b = <72/<7i coming 
from differences in the hopping amplitudes. Eq ^ implies 
a characteristic frequency scale 



Cai{d2 + bdi). 



(8) 



For reasons discussed earlier, in our discussion we will 
assume that the normal state conductivities cri_2(A = 0) 
are frequency independent and use these to define lu*. 

Fig 4 plots the Re and Im atuayer for the normal(A = 
0) state and for a = 0.1 superconducting state with 
LU* < 2A. The low-frequency suppression in the nor- 
mal state caused by the local field effects gives rise to 
a low-frequency upturn in r^^. The superconductivity- 
induced changes which have been discussed in detail for 
the single-layer case also apply to the bi-layer case. But 
new features may occur in the bi-layer case. In particu- 
lar, for some parameter ranges including those used in Fig 
4b, the local field effects lead to a bilayer plasmon (diver- 
gence in auiayer) at w < 2 A. This produces an additional 
zero crossings in Imcr, shown as the lower-frequency circle 
in Fig 4b. The second zero-crossing at a higher frequency 
(again circled) has it's origin in the divergence at the gap 
edge discussed earlier for the single-layer case (compare 
Figs 3c and 4b). 
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Fig. 4 Real and imaginary parts of tlie conductivity of a bi-layer 
superconductor with b = 0.075 calculated as described in ref |l2| 
for the normal state (A = 0) and for an a = 0.1 superconducting 
state with ui* = 0.4A. 

Corresponding to the two zero-crossings (Fig 4b) there 
are two delta-functions in Recr~^. These two delta- 
functions are obtained also for other values of a but their 
relative weights vary with the value of a. The form of 



normalized to 



oo limit for the case when the 



bilayer plasmon lies inside the gap for all non-zero values 
of a, coincides exactly with that for the corresponding 
single-layer case as shown in Fig 1. 
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Fig. 5 Optical scattering rate (Eq|l|) for a fei-Jaj/er super- 
conductor with b = 0.075 and an uj* = lOA normalized to — ► oo 
limit, in normal state (dashed line) and for superconducting state 
for different values of fluctuation parameters a = 0.1, 0.7 and 1. For 
superconducting state only the contribution for uj > 2A is shown. 

The detailed structure of zero-crossing and the form 
of conductivity varies with the ratio of the gap to the 
cross-over frequency lj* (Eq thus changing the form 
of T^^. For uj* >> 2 A, and a — 1 case for example, the 
bi-layer structure introduces a zero-crossing of Imcr just 



above the gap, which combined with the smallness of the 
Recr gives rise to a strong divergence in as shown in 
Fig 5. 

The sum-rule arguments can be easily extended to the 
bi-layer case and we get the sum-rule as given by |^ for 
the single-layer case but with f{a) replaced by 



f bilayer (o^) 



A7r(l 



di+bd2 2cti(A = 0) 



(9) 



In conclusion, we have studied the behavior of "op- 
tical scattering rate" Recr"^ in a model for the c-axis 
conductivity of high temperature superconductors. The 
presence of quantal and thermal fluctuations in the phase 
of the superconducting order parameter is shown to lead 
to a non-intuitive increase in the integrated area of the 
"scattering rate", Rea~^(u;). This increase is a net re- 
sult of the contributions from the lu < 2A and ui > 2A 
regions which depend sensitively on the strength of the 
phase fluctuations. In the strong phase fluctuation limit, 
Ke(T^^{uj > 2A) is decreased below its A = value 
whereas for weak phase fluctuations it is increased. While 
the variations in Recr"^ are seen to be a useful probe of 
order parameter phase fluctuations, the large contribu- 
tions at w < 2A coming from zero-crossings of Ima sug- 
gest that caution is required in interpreting Recr""'^ as a 
scattering rate. 
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